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CREEP-BUCKLING ANALYSIS OF RECTANGULAR-SECTION COLUMNS 


By Charles Libove 
SUMMARY 


A previous analysis of the creep behavior of a slightly curved pin- 
ended H-section colum under constant load is extended to the slightly 
curved solid rectangular-section column. The analysis leads to a dif- 
ferential equation for the plastic strains at the midheight cross 
section. The form of the equation indicates the significant parameters 
which may be useful in plotting test data on the creep life of columns. 
These are a lifetime parameter t'er, an initial-straightness param- 


eter S or S', and the ratio of the average applied stress to the Euler 
stress o/og. A numerical method of solving the differential equation, 
suitable for use with a high-speed digital computer, is described, and 
typical computed results are given. The existence of a finite lifetime, 
although not evident from the differential equation, is argued intuitively 
and confirmed by the numerical computations. 


INTRODUCTION 


The high temperatures that can develop in aircraft during supersonic 
Plight make it important to consider the possible Limitations due to 
creep on the useful service life of the structural components. Ina 
previous paper (ref. 1), the creep of a slightly curved pin-ended ideal- 
ized H~section column under a constant load and constant temperature was 
studied theoretically. The material of the colum was characterized by 
an assumed creep law, for constant uniaxial compressive stress and con- 
stant temperature, of the form 


€ = aa AeP OK 


tj 


where e is the total compressive strain, o is the constant compressive 
stress, t is the time after application of the stress, and E, A, B, 
and K are materiai constants. This form was selected because it 

applied to at least two alloys: the 75S-T6 aluminum alloy at 600° F 

(ref. 2) and a low-alloy steel at 800° F and, possibly, 1,100° F (ref. 3), 
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Shanley's engineering hypotheses of creep (ref. 4) were used in order to 
generalize the assumed constant-stress creep law to cover situations in 
which the stress varies with time, the condition encountered in the 
fibers of a column undergoing continuous lateral deflection because of 


creep. 


In the present paper, the analysis of reference 1 is extended to 
the solid rectangular-section column. Except for the shape of the cross 
section, the assumptions of reference 1 are retained. In addition to 
those cited previously concerning material behavior, these assumptions 
are: (a) that the initial shape is half a sine wave, (b) that the load 
on the column is applied rapidly enough so that negligible creep occurs 
during the loading period (that is, elastic behavior on loading) but 
not so rapidly that dynamic effects have to be considered, (c) that the 
shape of the deforming column remains at all times sinusoidal, (d) that 
plane sections perpendicular to the (curved) center line of the column 
before loading remain plane and perpendicular to the center line after 
loading, and (e) that the lateral deflections are small compared to the 
length of the column. 


The present analysis is very similar in its basic assumptions to 
that of Higgins (ref. 5) which also deals with the solid rectangular- 
section column; however, the two developments are entirely different, 
with the present development making evident the significant parameters. 


SYMBOLS 
A, B, B, K material constants 
b thickness of column 
Ij, I2 certain integrals, across column PRG of functions 
of pu at midheight 
L length of column 
n number of equal spaces into which cross-sectional thick- 


ness is assumed divided 


P | load parameter, en ae 


ef fal 


Ry functions of the values of u,, defined by equation (22) 
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ter 


1 
ter 


Uys Hos » 


Q] 


straightness parameter, e S 


alternate straightness parameter, œe 
time 


colum lifetime, value of t at which lateral deflec- 
tion become infinite 


1/K aa 
time parameter, (=) / eBo/K 
1/K 
AEB / Bo/K 
lifetime parameter, tar = e 


width of column 


distance measured from center line of column, positive 
toward convex side 


coordinates z of centers of the n spaces (starting 
from concave side) 


amplitude of initial sinusoidal deflection 
midheight lateral deflection (not including &ọ) 
compressive strain 

average compressive strain at midheight of column 


plastic-strain parameter, ral a 3 


values of u at centers of spaces 

dummy coordinate representing z 

compressive stress 

average compressive stress on column 
2 


x E 


Euler buckling stress of colum, ROEE 
12(L/b) 
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ANALYSIS 
Assumed Law of Material Behavior 
Because the longitudinal stress in any fiber of a continuously 
deflecting column varies with time, the description of the material 


given by the assumed constent-stress creep law mentioned in the 
introduction 


a = g + Ae BOE (ir) 


t 


is insufficient for purposes of analyzing the creep behavior of a 
column. Shanley's engineering hypotheses of creep (ref. 4) can be used 
to derive the following generalization of equation (Lis which covers 
behavior under varying compressive stress: (See pp. 460-461 of ref. 1 
for the derivation.) 


Ms 
H 
ki] Qe 


= f(o,€) | (2) 


where PA 


1/K 
e(aje) = EUe) (3) 


and the dots denote differentiation with respect to time. 


Equation (2) implies that the material behaves elastically unđer 
an infinitely rapid increase or reduction in stress. If equation (2) is 
rewritten with the term ô/E transposed to the left-hand side, f(c,¢) 
is seen to represent the creep rate, or rate of growth, of the plastic 


strain e - J, 
E 


Equation (3) implies that, if after a period of creep under com- 
pressive stress the stress is reduced to tension (a negative), the 
creep rate will still remain positive - that is, compressive - although 
it may become exceedingly small as the stress takes on larger tensile 
values. The neglect of tensile creep will introduce some unconservatism 
into the analysis when the lateral deflection gets so large that the 
stress becomes tensile on the convex side of the column. The degree of 
unconservatism, however, should be small in many practical situations 
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in which the initial curvature is reasonably small and the average 
applied stress is reasonably high but not close to the Euler stress. 

In such cases the time between the occurrence of appreciable tension 
and the calculated ultimate collapse of the columm will be a small part 
of the total column lifetime. 


Qualitative Description of Column Behavior 


The column, before application of load, is shown in figure I(aj 
The cross section is rectangular, of thickness b and width w. A 
slight initial curvature is present in the form of half a sine wave of 
amplitude &ọ. 


A load, less than the Euler load, is instantaneously applied and 
produces an average stress o and an instantaneous static deflection 
8,(0), as shown in figure 1(b). By hypothesis, the material behaves 


elastically during instantaneous stressing; therefore, the conditions 
existing immediately after load application (t = 0) are obtainable from 
an elastic analysis. At that instant no plastic strain exists anywhere 
in the column. 


With the load held constant, the column will continue to deform 
because of creep and will deflect laterally because of larger creep 
strains on the more highly stressed, or concave, side than on the convex 
side (fig. 1(c)). The problem is to determine the history of the 
deformations. 


That the present analysis should predict a finite lifetime for 
nonvanishing initial curvature can be established by the following 
intuitive argument: A solid rectangular-section column is certainly 
less stiff than the H-section column formed from it by concentrating 
half of the material of the column at each extreme fiber, all other 
conditions remaining unchanged. According to the analysis of reference 1, 
however, which is based on the same assumptions as the present analysis, 
the lateral deflection of any slightly curved H-section column approaches 
infinity in a finite time; hence, the deflection of the less stiff solid- 
section colum must certainly become infinite in finite time. 


Development of Basic Equations 


For a rigorous analysis, equilibrium, compatibility, and the law 
of material behavior must be satisfied at all sections of the column at 
every instant. For simplicity, however, the deflection shape is assumed 
to be merely a magnification of the initial shape - that is, half a sine 
wave. Because of this reduction to a single degree of freedom in the 
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deflection shape, the requirements of equilibrium, and so forth, can at 
best be satisfied only in some average way over the length of the colum. 
An even simpler solution, that followed herein, is to satisfy the require- 
ments at only one section of the colum, for example, at the midheight. 


Figure 2 shows the upper half of the column of figure 1(c) as a free 
body. The compressive stress o at the cut section varies from fiber 
to fiber and also varies in time - that is, o = o(z,t). Similarly, the 
compressive strain e = e(z,t). 


Equilibrium of vertical forces requires that. 


b/2 _ 
i o dz = ob (44) 
-b/2 : 
Equilibrium of moments requires that 
bp /2 
J oz dz = -db(& + 51) (5) 
-b/2. 


The assumed law of material behavior (eqs. (2) and (3)) gives the fol- 
_Lowing equation for the strain rate: 


1/K 
gd, Kaeo) (6) 


Finally, the assumptions that the shape is sinusoidal, that plane sec- 
tions remain plane and perpendicular to the center line of the colum, 
and that the deflections are small result in the following relationship 
between the strains and the lateral deflection at the column midheight: 


2 
x-6] 
Le 


E€ = -Z zee 


L a om, e s 
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where € is the average strain across the section, or 
€ = (€) 5-0 (8) 
and op is the Euler buckling stress of the colum; that is, 


nPE 
oo mae a (9) 
i 12(L/b)* 


Reduction of Basic Equations 


Equations (4) to (7) can be reduced to a single equation in terms 
of the elastic strain € - = at the midheight of the column, which is 


also a function of position and time. Equations (4), (5), and (6) ar 
first rewritten as follows (with cognizance being taken of eq. (7)): 


b/2 7 
| (e - gjs- Ps € dz 
-b/2 E B  J-b/2 


= b/2 120,83 m 
EE a e - a n 
= e(z - z) (10) 
M (< - Sz dz = Gae : ea) im Ez dz 
~b/2 E -b/2 
san + d; b/2 =- an 5 
a d es z dz 


H 


obs (0 - G\bs | 
O E 1 
T (11) 


E 
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1/K.Ba/ 
-E | aman 


ga o\K~* 
(6) 
120625 
=a -(e- g)te-2_ et 
1/K ap 
KA’ e 
- — (22) 
o\K 
E ee 
(6° § 


Equation (12) can be reduced further if € is eliminated through.use 
of equation (10) and 8 eliminated through use of equation (11). The 
result, after slight simplification, is the following equation in which 


€ - 5 is the only unknown: 


E aer oat gal (dee 
b/2 
T L g Fa. 3) wa ( j =)" “ A 
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where € is a dummy variable representing z. Introducing the following 
parameters; 


the plastic-strain parameter 


EB g\ _ 
Ze - g) =u ol 
the time parameter 
L/K 


the initial-straightness parameter 


_68¢ °0 IE 
K b gio 
e E = g (16) 
and the load parameter 
t 
=P (17) 
12 8 
OF 
permits equation (13) to be simplified to 
1/2 1/2 
(| ; Eal 2 4$) 
D ilk. gezh,\ ei? ated (18) 


ot? 


where p =u(z,t) is subject to the initial condition u(z,0) = 0. 


Equation (18) gives, essentially, the rate of growth of the plastic 
strain to the 1/K power for any fiber of the colum as a function of 
the position z of that fiber, the instantaneous value of the plastic- 
strain parameter u in that fiber, and the two integrals involving the 
instantaneous values. of u in all the fibers of the colum. 


LO NACA TN 2956 


For the initially perfectly straight column (S = 1), equation (18) 


Will be satisfied by yu = tK, as is expected, inasmuch as u = t'K isg 
equivalent to equation (1), the law of material behavior for direct 


compression. 


Numerical Method of Solution of Equation (18) 


Equation (18) lends itself readily to numerical solution by means 
of the modified Euler method (ref. 6). The following adaptation of this 
method has been found suitable for use in conjunction with a high-speed 
digital computer, the National Bureau of Standards Eastern Automatic 
Computer (SEAC): 


Let the cross section of the column at the midheight be divided into 
m equal spaces. The width of each space will be b/n and the centers 
of the spaces, starting at the left-hand, or concave, side of the column 


and going toward the right, will be located at z = 2, = - 5 + s b, 
n 
Z = 25 =- 24 2R, . «© «. Let the values of the plastic-strain param- 
n 


eter p at the centers of these spaces be denoted by Hi; Ho, * * * Mp? 


respectively. The two integrals appearing in equation (18) can then be 
approximated in terms of the parameters wu; and Z,, the very simplest 
approximation for machine-calculation purposes being that obtained by 
assuming that the integrand is constant over the width of a space and 
equal to its value at the center of the space. Then, 


| 


1/2 
Lp t ESTEN aterse un) (19) 


a 
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Writing equation (18) for the specific z = Zj; Zo; * + + Zņ gives 


the following system of ordinary differential equations which is approxi- 
mately equivalent to equation (18): 


it 
al Ly G pr 1 -py+Ty +6 (-1 + Ler, 
= e 


dt' 


a uE E 2 | uptty6 (-1+ 22A 
aoe (21) 


1/K B en-1 
a(i, / ) a 3 [nun +6 (-24 = BA 
dt! p 
Letting 
eee a pussy +6 (-1+ 21-1) Pt] (22) 
Ry = § e 


simplifies equations (21) to: 


a u L/K 

att o a 

1/K 
ove) = Ry (23) 
a(u,2/*) ie 


le NACA TN 2956 


The parameters Ry are functions of t' because they involve the 
ps (t'). 
Equations (23) define curves of the py plotted against t' and 


can be solved in a step-by-step fashion subject to the initial conditions 
u4(O) = 0. As may be expected on the basis of the physical argument 


presented previously, at least one of the pui, probably wy, will tend to 
approach infinity asymptotically at some finite value of t'. A step- 
by~step solution of equations (23) on the basis of equal increments 

of t' therefore should not be attempted. Rather, the solutions should 
be carried forward on the basis of equal increments in p] or u LK, 


and the increment in t' should be regarded as one of the unknowns in 
each step. A detailed account of the calculation procedure follows. 


Suppose that the curves of His Hos * * * Hy plotted against t' 
have been traced up to some particular value of t'. It is now desired 
to carry these curves forward one more step corresponding to a specified 
finite increment in pu} denoted by Ay . To be determined are the 
increment At' in t' and the increments Ao, Aiz, + + + An in the 
other uy corresponding to the specified increment Au, in yy. Since 


the values wu, and Ai, are known, the increment A(u,2/K) in the 


quantity „1/5 can be computed. The first of equations (23) then 
gives the following approximation to At': 


t 1/K\ 1 
a = Afat FEOT (24) 


With At' known, the remaining equations give the increments A(ug!/K), 
A(us/*), Eo ie <8 afun 1/£) as follows: 


Afua?) = atnl) 
a(us1/5) = At'Rs(t") 
, (25) 
A(u,M/K) = at'R, (t") 
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From these increments, the values of po, Has ie Hy corresponding 


to t' + At’ can be computed. The curves have now been tentatively 
advanced one step and tentative values of R,(t' + At') can be deter- 


mined. Improved values of At! and alu 1K), A soe o alunt) 
can now be obtained by replacing the right-hand sides of equations (2}4) 
and (25) by averages based on the values of R; at the end of the 


interval as well as at the beginning of the interval, that is, by using 
the equations 


afu t/) = At! Lla) + Ro(t' + ao) 
Auz) = At' z s(t") + Ra(t' + ae") 
(27) 


alt) = At! z}aa(t") er cue at) 


The improved values af u,(t' + At') obtained from equations (27) 


give new values of R,(t' + At') which can be substituted in equations (26) 
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and (27) in order to obtain still better values of At' and A(up1/K), 


Ae so. Alun t5). This process can be repeated until no change 


is evident in the values of At' and u,(t' + At’), palt" + Ab'), ee 
u(t! + At’). 


The procedure just described can also be used without modification 
to start the curves. It is necessary to use only the fact that wu y(0) = 0. 


Approximate Solution of Equation (18) Based on n = 2 


If, for purposes of numerical integration, the cross section Were 
divided into two equal spaces only, then the system of equations (21) 
would reduce to the following pair of equations: 


a(u,1/K _ -1/2 (#172) ee) 
dt” 


(28) 


With the differentiations indicated on the left-hand sides carried out, 
these equations become identical in form with the equations which arise 
in the analysis of the two-flange column - namely, equations (21) and (22) 
of reference 1. The solution presented in reference 1 can therefore be 
used to obtain the approximate solution of the present problem corre- 
sponding to n = 2., In order to use the solution in reference l, certain 
substitutions should be made in the symbols and graph labels of that 
paper and, then, the results of that paper will apply directly to the 
present problem. The substitutions are as follows: 


(1) Replace Yp/71, by 6 


(2) Replace ny by un (FE - z) 
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(3) Replace ng by u4 > =) 


L T 1 
(4) Replace 7,t by (= - i) [5s -1/24 


P 1/K 
(5) Replace Ypter by (= z i) Ps 1/24 aT 


(6) The constant K is unchanged 


Lateral Deflections, Strains, and Stresses 


With the yu 4(t') completely determined by the solution of equa- 
tion (18), the time history of the plastic strains e - = is, in effect, 
known. Equation (11) can then be used to compute the midheight lateral 
deflections 6& ‘coe and equation (10) gives the history of the average 
midheight strain e(t). The strains e(z,t) are then obtainable fram 


equation (7). A knowledge of e - = and e is sufficient to determine 
the midheight stresses o(z,t). 


RESULTS AND DISCUSSION 


Significant Parameters 


The usefulness of data obtained from tests is enhanced if the data 
are presented in terms of the appropriate parameters. Equation (18) of 
the present analysis suggests that test data on the creep life of coluwms 
whose material creep curves are adequately represented by equation (1) 
may profitably be plotted in terms of the parameters K, tl op? p 
and @/cp;, Where t',, is the value of t' corresponding to collapse of 
the column. If the analysis were perfectly applicable, the test data 
would give a single curve of t'ar plotted against S for a given value 
of o/OR and a given value of K. Scatter of test points away Prom such 
a curve will, of course, occur, because of shortcomings of the theory - 
for example, the neglect of third-stage creep in equation (1), the approxi- 
mate engineering hypotheses used to generalize equation (1), the neglect 
of tensile creep in the generalization, the assumption of an exactly 
sinusoidal shape for both the initial and the subsequent deflections, 
and the assumption of instantaneous elastic loading. 
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Inasmuch as the parameter S contains the parameter C/o and each 
plotted curve or test point has a definite value of o/oR associated with 
it, a simpler parameter S' may be used instead of 5, where 

6Ba ĉo 


a1 ., kK D 


z 
-s 3 (29) 


If the test results for a single material at a single temperature 
are to be plotted, then A, B, E, and K are the same for all the 
tests and the significant parameters may be simplified further to 


Bo 5 bs 
Tags o/K oo and 0 / Op 


Calculated Values of Dar 


In the section entitled "Analysis" a numerical method of solving 
equation (13) was given, based on dividing the column cross section at 
the midheight into n equal spaces and writing equation (18) for the 
center of each space. This method has been carried out and found to be 
satisfactory for the case K= 1l, 0/oR = 0.7 and 0.9, and for several 
values of S. The value n = 1} was found to lead to sufficiently 
accurate results in a trial calculation, and this value was used for 
all the cases. 


Figure 3 shows a typical curve of p} plotted against t'. (The 
curves of Ho; Uz, * * * Hy are not shown, but they would lie suc- 
cessively below the curve of u1.) This figure tends to confirm the 


intuitive argument presented previously to the effect that the deflec- 
tions should approach infinity in a finite time. Computation was stopped 
when the computed increment in t' corresponding to the specified 
increment in puq became less than the precision of the computation. 


The last computed value of t' was taken as being equal to the asymptotic 
value t'er for all practical purposes. E 


The values of t'or taken from a number of curves like the one in 
figure 4 are summarized in figure 4 in plots of t'er against S for 
a/e = 0.7 and 0.9. The circles are the actually computed data through 


which the curves are faired. 


As is to be expected, t'er tends toward zero for columns that are 


more and more initially curved (S—~ 0) and toward infinity for columns 
approaching straightness (S—»>1). In order to show more clearly the 
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£ behavior of the curves in these two regions, the data of figure 4 are 
replotted in figures 5 and 6 in which logarithmic scales are used to 
magnify the region of S = 0O and S = 1, respectively. 


All the computations were performed on the National Bureau of 
Standards Eastern Automatic Computer. The computation of a value of 
t'er took about 15 minutes, exclusive of the exploratory computation 
time required to establish optimum values of precision, the step size 
A; and the number of cross-sectional spaces n. The calculated values 


of t'er are felt to be accurate to within 2 percent. 


An approximate solution for the solid rectangular-section column 
corresponding to n = 2, as is indicated in the section entitled 
"Analysis," already exists implicitly in the two-flange-colwm analysis 
in reference 1. Figure 5 of reference 1, which summarizes the theoret- 
ical lifetime data for two-flange columns, can therefore be used to yield 
approximate information on the lifetimes of solid columns by replacing 
the label IRITI by the label S and the label Yrter by the label 


-1/2 
- (F zy Ss ot 
4 2 


er 


In figure 7, the approximate results thus obtained for the 
rectangular-section column for the cases where K = L and o/oR = Ouse 
and 0.9 are compared with the much more accurate results for these cases 
based on computations with n = I} and obtained from figure 4 of the 
present paper. 


CONCLUDING REMARKS 


On the basis of certain assumptions regarding the shape of a column 
and the creep behavior of the material under varying stress, an analysis 
has been made of a slightly curved solid rectangular-section pin-ended 
column carrying a’constant load. The analysis led to a differential 
equation for the plastic strains at the midheight cross section as a 
function of time. The form of the equation indicates the significant 
parameters of the problem, which parameters may be useful in plotting 
test data on the creep life of columns. These are a lifetime parameter 
t'er; an initial-straightness parameter S or S', and the ratio of the 
average applied stress to the Euler stress G/o,. 


A numerical method of solution of the differential equation, suitable 
ee — for use with a high-speed digital computer, has been described, and. 
typical computed results have been presented. 
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The existence of a finite lifetime, although not evident from the 
differential equation, was argued intuitively and confirmed by the 


numerical computations. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., March 31, 1953. 
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Figure 1.- Schematic description of colum. NACA 
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Figure 2.- Upper half of column of figure l(c). 
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Figure 3.- Variation of plastic-strain parameter Hy, with time 
parameter t'. K = 1; c/o, = 0.7; 8 = 0.1. 
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Figure 4.- Variation of lifetime parameter t'er With initial-straightness 
parameter S for two values of the load parameter o/op end Tor Ked; 
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S 


Figure 5.- Replot of figure 4 showing more clearly the nature of the curves 
in the region of S = 0. 


NACA TN 2956 23 


Figure 6.- Replot of figure } showing more clearly the nature of the curves 
in the region of S= 1l. 
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Figure 7.- Comparison of approximate solution, based on n = 2, with 
accurate solution, based on n = l4, for K= 1. 
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